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Abstract. Work in the paradigm of quantum fluctuation theorems of Crooks and 
Jarzynski, is determined by projective measurements of energy at the beginning and 
end of the force protocol. In analogy to classical systems, we consider an alternate 
definition of work given by the integral of the supplied power determined by integrating 
up the results of repeated measurements of the instantaneous power during the force 
protocol. We observe that such a definition of work, in spite of taking account of 
the process dependence, has different possible values and statistics from the work 
determined by the conventional two energy measurement approach (TEMA). In the 
limit of many projective measurements of power, the system’s dynamics is frozen in 
the power measurement basis due to the quantum Zeno effect leading to statistics 
only trivially dependent on the force protocol. In general the Jarzynski relation is not 
satisfied except for the case when the instantaneous power operator commutes with 
the total Hamiltonian at all times. We also consider properties of the joint statistics 
of power-based definition of work and TEMA work in protocols where both values 
are determined. This allows us to quantify their correlations. Relaxing the projective 
measurement condition, weak continuous measurements of power are considered within 
the stochastic master equation formalism. Even in this scenario the power-based work 
statistics is in general not able to reproduce qualitative features of the TEMA work 
statistics. 
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Transient fluctuation theorems are exact relations restricting the statistics of work 
performed by externally controlled classical forces. While the considered system initially 
must be in a thermal equilibrium state, the subsequent forcing may drive it to out 
of equilibrium to states that cannot be described in terms of linear response theory. 
Exact relations of such form were pioneered by Bochkov and Kuzovlev [I] and the ones 
pertinent to our discussion bear the names of Jarzynski |2| and Crooks [31. The Crooks 
relation [3|, 

Pa(w) = e~P {AF ~ w) pA(~w) (1) 

relates the probability density function (pdf), Pa(w), of work performed on a system 
with Hamiltonian 77[A(f)] by the action of a generalized force \(t) that varies according 
to the protocol A = {A(t)|0 < t < r} to the work pdf of another process governed 
by the time-reversed protocol A = {A(r — t)\0 < t < r|§} Both processes start in a 
canonical equilibrium state at the same inverse temperature fl described by the density 
matrices p t = Z~ l (t)e~P H ^ x<yt ^ with t = 0 and r for the forward and the backward 
process, respectively. Here Z{t) = Tr for t — 0 and r yields the partition 

functions which determine the free energy difference A F = —fl~ l ln(Z(r)/Z(0)). 

A straightforward integration of (|Tj) with e~^ w brought to its left hand side leads 
to the Jarzynski equality |2jj, 

(e~ pw ) = e~ pAF . (2) 


Fluctuation theorems have been shown to be valid in a variety of situations ranging from 
open classical systems [4], closed j^-7| and open quantum systems |8—TO|. For recent 
reviews of this topic 11 121 (focus on classical systems) and 13 -16| (focus on quantum 
systems) can be consulted. 

These theorems have also been verified experimentally for a variety of classical 
systems IT 201 and have also been used to determine free energy changes in molecular 


systems by measuring work statistics in controlled non-equilibrium processes 17,18 


Further details on this topic are available in the reviews 21,22 


The situation regarding experimental verification of fluctuation theorems for 
quantum systems is still in its nascent stages. One central issue in the quantum context 
is that work is not an observable (7j; it is process dependent and therefore the potential 
values that it may take in a single experiment cannot be represented as the eigenvalues 
of a hermitian operator acting on the Hilbert space of the considered system. This 
issue underpins the most common approach to determine work statistics for a quantum 
system, the so-called two-energy-measurement approach (TEMA) [7,14,15 : here the 
work supplied by a time dependent force \{t) during the time interval 0 < t < r is 
defined as the difference between the system’s energy at the final and initial times. The 
TEMA can be divided into the following steps: (i) the system is prepared in a thermal 


§ For the sake of simplicity we restrict ourselves to forces A(f) which transform evenly under time- 
reversal. 
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equilibrium state at inverse temperature ft at t — 0 with the initial Hamiltonian i/[A(0)]; 
a projective measurement of energy is performed yielding one of the eigenvalues 


n 


e n (0) of //[A(0)] as possible outcome; due to the measurement the state of the system 
is projected to the associated energy eigenstate |n;0); (iii) subsequently it undergoes a 
unitary time evolution during 0 < t < r governed by the Hamiltonian H[\(t)] which 
changes in time according to the prescribed protocol A; (iv) at t — r a second projective 
measurement of energy is performed; its outcome is an eigenvalue e TO (r) of the final 
Hamiltonian H[\(t)]. The work done in such a process is given by w e = e m (r) — e n (0). 
In general, projective measurements turn out to be difficult to perform in an experiment. 
Nonetheless proposals for implementing such a scheme using trapped ions was suggested 
and its experimental realisation was recently reported in 24 . In an indirect 


m 


23 


verification of the fluctuation theorems detailed in 25^27 , the experimental difficulties 
of projective measurements were circumvented by encoding the characteristic function 
of work in the reduced state of an ancillary system that interacts with the system of 
interest with a strength determined by the force protocol. An experiment, using NMR 

. A second way to address the 


spectroscopy, of such a proposal was reported in |2 
difficulty of performing projective measurements is to ask if the replacement of the 
same by non-projective generalized measurements, which may be easier to implement in 


an experiment, preserves the fluctuation theorems 29 301 (see also (31,32 for a recent 
proposal to measure work as the outcome of a single generalised measurement). With 
this line of enquiry we found in an earlier publication 29 a no-go theorem restricting 


projective measurements of energy as the only ones within a large class of generalized 
measurements that satisfy both Crooks and Jarzynski equalities for arbitrary protocols. 
Although, for some specific forms of the generalized measurements, modified fluctuation 
theorems may be derived 


30 


In this work, we focus on yet another way to address the issue of measurement of 
work in quantum systems. Conventionally, in classical systems work can be determined 
in an incremental way by integrating up the supplied power which can be inferred from 
the instantaneous state of the system 21 . A direct extension of such a method to 


quantum systems is difficult since the system will have to be continuously monitored 
33,34 . It is known that a continuous monitoring of a quantum system can cause the 


freezing of coherent dynamics, a phenomenon known as the quantum Zeno effect 33,34 


The central aim of this work is to explore in detail how and why the extension of power 
measurements to estimate the work done on quantum systems fails to provide work 
statistics that satisfies the transient fluctuation theorems of Crooks |3| and Jarzynski (2|. 

The paper is organized as follows. In section [2] we define the quantum mechanical 
variable corresponding to the instantaneous power supplied to a system that experiences 
a force via a coupling to a generalized coordinate. We introduce an estimate of work 
computed by incremental projective measurements of the power in conjunction with 
the protocol, and explain how such an estimate differs from the usual TEMA based 
work both in magnitude and range of the possible values. In the limit of frequent 
measurements, the system’s unitary evolution is frozen due to the quantum Zeno effect 
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(we refer to such a regime of the dynamics as the “Zeno limit” henceforth). We analyze 
the statistics of the power-based work in this limit and point out how the commutator 
of the power operator with the total Hamiltonian is crucial in deciding the merits of the 
work estimate. We also derive an inequality for the dissipated work estimated by the 
power measurement in the Zeno limit. Further, the results so far obtained are illustrated 
by the example of a two-level system undergoing an avoided crossing. In section [3j 
we allow for energy measurements at the beginning and the end of the force protocol 
additionally to the coordinate measurements necessary for the power measurements 
and compare the distributions of TEMA, w e , and power-based work estimates w p . In 
section[4]we relax the assumption of projective instantaneous measurements of the power 
operator and adopt a weak, continuous measurement of power to define the work. The 
treatment there is carried out within the framework of Stochastic Master Equations 
(SME) [35,361. Work statistics obtained from numerical solution of SMEs for the two- 
level system example are also studied. In section [5] we conclude the paper. Appendices 
[A] and [B] provide additional details omitted from the main text of the paper. 


2. Power measurement and the quantum Zeno effect 


In classical mechanics, the energy-work relation can be invoked to define the work 
supplied to a thermally closed system[jj]as the increase in its energy. First, let us consider 
a classical system whose energy is determined by the Hamiltonian H cl [ z, A(f)] where z is 
a point in phase-space and A (t) is an external parameter which is varied in time leading 
to a change in the system’s energy. Hence the work done by changing A (t) in the interval 
0 < t < r according to some specific protocol is given by |2 

w cl = H cl [Z(r, z), A(t)] - H d [z, A(0)], (3) 


where Z (t, z) is the solution of Hamilton’s equation of motion 
|z = {H d [Z, A(«)], Z} 

at time t for a trajectory with the initial condition Z(0,z) = z. (J3]) can be rewritten as 
an integral of the total time-derivative of the Hamiltonian which agrees with the partial 
time-derivative. This gives an equivalent expression for work, 


w cl = I dt 


9ff d [Z(t,z),A(t)] 

dt 


where the integrand is the instantaneous power L cl supplied to the system at time t, 


L d (t, z) = 


[Z(t,z),A(t)[ 

dt 

OH [Z(f, z), A(f)] d\(t) 


d\ 


dt 


( 4 ) 

i.e. 

( 5 ) 

( 6 ) 


| We denote a system as thermally closed if it does not exchange energy and/or particles with its 
environment. In order to be able to perform work on it, it must though be possible to externally 
change some of the system parameters. 
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The power-based work expression (|4]) is more convenient to determine work in 
experiments |l7j. Since the initial condition z is only specified in terms of a probability 
distribution (for, e.g., a canonical distribution if the system is initially in thermal 
equilibrium), work becomes a random quantity. 

In quantum mechanics, both expressions ([3]) and Q can in principle be extended 
to operational definitions of work which, however, turn out to be no longer equivalent 
to each other. Here, we note that there have been earlier attempts 37-42 to define 
work in quantum systems using expressions analogous to (J4]) . In 137 this is done in the 
context of a driven harmonic oscillator interrupted by a small number of measurements 
of the coordinate and in 40 the same system is examined under the continuous quantum 
histories framework for the power operator. In the context of superconducting Cooper- 
pair box systems some lower order moments of a power-based work were considered 

In some 


m 


41 for a closed system and extended to include open systems in 42 


of these attempts 38,41], power-based work has been treated as a standard quantum 
mechanical observable (with a corresponding operator) in contrast to the view point 
in |7|, which we adopt. One central aim of this work is to show that even a more careful 
implementation of the power-based work for quantum systems taking into account 
the process dependence will generally lead to qualitatively different statistics from the 
TEMA definition of work and generally fail to satisfy the fluctuation relations of Crooks 
and Jarzynski. 

To that end, we consider the simplest possible situation where a single scalar 
parameter A (t) acts as a force on a system via a generalized coordinate A". In the 
absence of this force the system is described by the Hamiltonian Hq. In the quantum 
scenario, the generalised coordinate X corresponds to a self-adjoint operator X acting 
on the system’s Hilbert space. The total Hamiltonian is then given by: 

H[X(t)} =H 0 + A (t)X. (7) 

In analogy to the classical form of power ([b]) we define the power operator L(t) as 


m = 


dH H [\(t)\ 


dt 




( 8 ) 


where the superscript H denotes the Heisenberg picture and the dot a time derivative. 
Since the generalised coordinate does not explicitly depend on time, we simply indicate 
the Heisenberg picture by the presence of the time-argument omitting the superscript H. 
In order to determine the power, one needs to perform a measurement of the generalized 
coordinate A". The possible measurement outcomes, in a projective measurement 43 , 
are determined by the eigenvalues x a of the generalized coordinate X = Yh a x °^a- For 
the sake of simplicity we assume that X has a non-degenerate discrete spectrum. Hence, 
the eigenprojection operators = \<p a ) (<p a \ are determined by the eigenfunction \ip a ) 
of X. In order to precisely capture the work as the integral of the power one should, in 
principle, continuously measure the generalized coordinate. This though will inevitably 
freeze the dynamics of the system in an eigenstate of X. Provided that the force protocol 
is sufficiently slow compared to the unitary dynamics of the system, one can try to avoid 
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the full halt of the dynamics by performing only a finite number of power measurements 
and approximate the integral by a discrete sum: 

N 

w p N) = ^2 X(ti)x ai h, (9) 

2=1 

where we assumed that N measurements take place at regularly spaced times t, with 
ti- |_i — ti = t\ = h in the interval (to = 0) < t < (tv+i = t). The first basic difference 
between the power-based and TEMA work estimates emerges from the varying allowed 
values of work that the two approaches produce. On the one hand, the set of possible 
TEMA work values is given by 

W e = {w = e m (r) - e n (0)|m G /(r),n G J(0)} , (10) 

where I(t ) is the set of indices labelling the spectrum of H(t ) = Ylneilt) e n(t)H n (t)- Here 
U n (t) are operators projecting on the eigenstates corresponding to e n (t). Comparing 
with all possible power-based work values having the form given by (J9]) , it is apparent 
that both the number of possible values and magnitudes of work are different in the two 
approaches. Indeed for a system with finite -D-dimensional Hilbert space, the maximum 
number of possible work values from TEMA is D> 2 , whereas in the power measurement 
case with N measurements it is given by D N . Whereas the allowed TEMA work values 
are functions of all parameters entering the full Hamiltonian ([7]), the power-based work 
values depend solely on the time-derivatives of the force at the measurement times. Only 
the probabilities which characterize the occurrence of these work values may depend on 
the other parameters of the system. 

Still one could hope that the two estimates of work might be similar in a statistical 
sense. For this purpose we now compute the distribution for the work estimate 
in 0 to compare it with the distribution of the TEMA work values in terms of 
their cumulative probabilities. Beginning at t — 0 with a canonical density matrix 
po = Z ~ 1 (we obtain for the joint probability £P\(~x) of finding the eigenvalues 
x = ( x aii x a 2 i • • •, x a N ) hi the N measurements at times 4, t 2 , ■ ■ ■ An- 

^a(x) = Tr 14 (x) poV\ (x), (11) 

where 


14 (x) = nyivny t • ■ • v&Zjh. (12) 

Here 14 = U tkt t k _ 1 ( A) denotes the time-evolution operator for the system from time 4 -1 
to 4 under the influence of the protocol A. It is the solution to the Schrodinger equation 
r)U 

ih= H[\(t)}U t , s (A), (13) 

obeying the initial condition 


EM A) = 1. 


As evident from (12), we have taken the measurements of the generalized coordinate 
to be projective (following von Neumann |43|). We will relax this assumption in 
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Appendix [A] and consider generalised measurements of the coordinate. Expressing the 
projection operators in terms of the eigenfunctions of X, the joint probability becomes: 

TV—1 

^A(x) = Y[ \ {Va k+1 \Uk+l\<Pa h )\ 2 (<P ai \UipoUl\(p ai ) . (14) 


k =1 


The pdf p^J (w) of finding the value w for the power-based work estimate Wp N ^ is then 
given by: 

N 





x ai h 


(x). 


(15) 


1=1 


For a large number of measurements, the Zeno effect 33,34 causes a freezing of the 
system in the state corresponding to the outcome of the first measured eigenvalue of A". 
This is a consequence of the fact that the transition probabilities |(<£ afc+1 | E4+i |<A* fc )| 
between different eigenstates of X vanish as h 2 . Hence for large values of N (which 
correspond to small values of h = r/(N + 1)), the joint probability in (14) becomes 


TV-1 


^a(x) = JJ S x 


*fc + l 


(wlpo I<A*i) + 0(h). 


k =i 


Putting this asymptotic result into the expression in (15) we obtain 
Pa”’W = Yi S - I A ( T ) “ A (°)] X «) W Pa Wc.) , 


(16) 


(17) 


for the work pdf neglecting corrections of the order h. Note that in the above h —> 0 
limit, the sum in the delta function specifying the work in (15) can be replaced by an 
integral yielding lim^oo YlZ=i KU)x ai h = [A(r) - A(0)] x ai . 

Having obtained the pdf for the power-based estimate of work, we can immediately 
check if the Jarzynski equality is satisfied. To that end we have for the average of the 
exponentiated work the following expression: 

r tv p -/3[a(t)-a(o)]x ~p\h 0 +\{o)x\ 

(e~^) r = I dwp%(w)e-e« = —- ^ -. (18) 

For the Jarzynski equality to hold, the numerator of the right-hand side must coincide 
with the partition function Tr . However, this is only the case if the unperturbed 

Hamiltonian H$ commutes with the generalized coordinate X, i.e. if [Hq,X\ = 0, which 
is an atypical situation. Because in the commuting case, the forcing does not lead 
to transitions between different eigenstates of H[X(t)], the set of allowed TEMA work 
values is given by 

W e = {(A(r) - A(0))x a |a G /} (19) 

and hence coincides with the set of allowed power-based work values. This follows from 
the eigenvalues of H [A(t)] taking the form e a (t) = e a + A (t)x a , a G / where / is the 
time-independent set labelling the eigenvalues e a and x a of H 0 and X, respectively. 
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In general, when Hq and X do not commute, the power-based work estimate does 
not conform with the Jarzynski equalitjj^J In this case one finds with the Golden- 
Thompson inequality 44 Tr e A e B > Tr e A+B (with A and B being hermitian) in (18) 
yielding 

Tr p -/3[A(r)-A(0)]X-/3[H 0 +A(0)A'l 

<e-*% > —---- 


m 


= e~^ F . 


Thus the power-based work estimate satisfies the inequality 

( e -^ w -AF)) p > i. (20) 

Hence, an estimate of the free energy based on the exponentiated power-based work 
average underestimates the true value 

AF p = -/T 1 ln(e-^) p < A F. (21) 


A more detailed comparison can be made on the basis of the power-based work pdf (]1T|) 

)df p Ae {w) given by [I 

e m (r) + e n (0 ))p A (m,n), 


and the corresponding TEMA work pdf p Ae (w) given by [14,15 
PAe(w) = Y' 5 (W~ 


( 22 ) 


where the joint probability p A (m, n ) to find the eigenstates |n; 0) and |m; r) in projective 
energy measurements at the beginning and the end of the force protocol, respectively, 
reads 

p A (m,n) = Tr U m (t) U Tja (A) U n (0)p 0 n n (0) U ] Tfi (A). (23) 

In the above equation n.,(f) denotes the eigenprojection operators corresponding to the 
eigenvalue eflt) of the Hamiltonian H [A(t)j and U T) q(A) is the unitary time evolution 


operator from the beginning to the end of the force protocol defined in (13). 


We further elucidate the differences of the work pdfs, (15) and (22), which both are 
rather involved, by means of a simple example. 


2.1. Landau-Zener 

Next we want to illustrate the differences between the power-based estimates and the 
TEMA work values and also the approach to the Zeno limit for a driven two-level system 
undergoing an avoided level-crossing. 

The Hamiltonian for a two-level system driven through an avoided crossing, known 


as the Landau-Zener(-Stiickelberg-Majorana) model 45 , is given by 


vt 

HlzW = Aa x + — (7 Z , 


(24) 


where Oj are the Pauli matrices. In this case the power operator is given by L — 
va z /2. Thus the power measurement involves projective measurements in the cr z basis. 
Digressing from our previous convention, we follow the usual custom and consider a 

In [40| a similar conclusion was obtained using the continuous quantum histories framework for the 
model system of a center-shifted harmonic oscillator. 
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symmetric protocol about t — 0 that runs between — r/2 < t < r/2. The Hamiltonian 
(24) is readily diagonalized yielding for the eigenvalues ej(t) = (—1) J y/ut/2) 2 + A 2 
with j = 1,2. Hence, the possible TEMA work values are given by: 


W e = {-E max ,0,E max }, (25) 

with E'max = 2 a/( ur/4) 2 + A 2 , whereas the set of possible work values based on N power 
measurements becomes 

WP= { 2 (JV+l) g ' » = - Ar .- iV + 2 ."'. iV }. (26) 

clearly showing a fundamental difference between the two approaches. We also note 
here that for this example the range of TEMA work bounds the possible values in the 
power measurement estimate. The range of the latter increases with the number N 
reaching ±ur/2 for N oc. The magnitude of the maximum work value possible from 
the TEMA, E max , is always larger than vr /2 but in the limit of large enough r such 
that vr /4 A, it approaches the latter. 

According to 0. the work pdf for the LZ problem estimated by power 
measurements in the Zeno limit becomes 


= ( z +\Po \ z +) K w ~ VT / 2 ) + ( z -\Po \ z - ) S(w + ur/2), (27) 

where \z±) denote the eigenstates of a z and p 0 = e^' !3HLz(y ~ T ^ T> /Z{—r/ 2) is the initial 
density matrix. On the other hand, the work pdf from the TEMA is given by |46| 

Pxeiw) = p e P LZ S(w + £ max ) + (1 - Plz)8(w) + PgPizbiw - E mSLX ), (28) 


denotes the LZ tunneling 


with p g = (l + e“ ,3Emax ) 1 , p e = 1 —p g and Plz = e~ 27rA2 P hv ^ 

probability [45] from the ground state at t — —r/2 to the excited state at t = r/2. 

Figures [T| and [2] display the cumulative probabilities of work computed for different 
numbers of power measurements N and for two different temperatures. The cumulative 
probability Q(w ) = Jfl dw'p{w ) quantifies the probability to find a work value 
w' < w based on the pdf p(w). As a quantitative measure between two work 
distributions we use the Kolmogorov-Smirnov (KS) distance D K $, corresponding to 
the maximum absolute difference between the respective cumulative probabilities, i.e. 
Dks = max,,; \Qflw) — Q 2 (if)|- For the parameter values chosen (see figure caption) the 
initial eigenstates of the LZ system have significant overlaps with the a z eigenstates, i.e., 


-r/2)) 


LZ 


| z±) (the minus sign on l.h.s stands for the ground state). As a result, in 


the low temperature example in figure [T| one can see that as N is increased, in agreement 
with (27), the largest jump of the cumulative probability occurs near w ~ ur/2. Only a 
small jump at w ~ —ur/2 is visible due to the low initial occupation of the excited state. 
For the high temperature example in figure [2j two jumps of comparable height appear 
in the large N limit as both eigenstates are occupied in the initial distribution. The 
cumulative probability from the TEMA is also plotted for reference and it is apparently 


+ Strictly speaking this expression is valid only for r —> 00 but provides a very good approximation for 
large finite r with vr/ 4> A. 
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Figure 1. Cumulative probability of work estimated from the sum of N power 
measurements for the LZ problem with v = 5A 2 /h and r = 20h/A (N = 10 - red, 
N = 10 2 - blue, N = 10 3 - green, and N = 10 4 - black with ascending line thickness). 
The initial temperature is small /3E max = 10. For large N , due to the Zeno effect, 
a distinct peak appears at the maximum value of w (see text for discussion). For 
comparison the pdf computed from TEMA is also shown (cyan dashed line). The inset 
shows the Kolmogorov-Smirnov distance between the TEMA and power-based work 
estimate as a function of N. 



Figure 2. Cumulative probability of work estimated from the sum of N power 
measurements for the LZ problem for the same system sweep parameters (and same 
legends) as in figure[l] The initial temperature is large /3E max = 10“ 4 . For large N, due 
to the Zeno effect, two distinct peaks appear at the maximum and minimum value of 
w (see text for discussion). For comparison the cumulative probability computed from 
TEMA is also shown (cyan dashed line). The inset shows the Kolmogorov-Smirnov 
distance between the TEMA and power-based work estimate as a function of N. 
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Figure 3. Cumulative probability of work from TEMA during a protocol for the LZ 
problem that is interrupted by N measurements of a z . The system sweep parameters 
and legends are as in figure[l] in particular, the color code indicates the same number of 
power measurements N. The initial temperature is small /3-E max = 10. For comparison 
the cumulative probability computed from N = 10 4 measurements of power is also 
shown (cyan dashed line). The inset shows the Kolmogorov-Smirnov distance between 
the TEMA and power-based work estimate as a function of N. 



w/A 

Figure 4. Cumulative probability of work from TEMA during a LZ sweep that is 
interrupted with power measurements. System parameters (and legend) are as in 
figure [3} The initial temperature is large j3E max = 0.1. 


quite different from the pdf for the power-based work estimate. The largest Kolmogorov- 
Smirnov distance between the TEMA and the power-based work distributions results 
in the Zeno limit. 


In earlier work 49 , it has been shown that though fluctuation theorems for TEMA 
work are robust to measurements during the protocol, the work statistics itself can 
be strongly modified. In this light it is interesting to compare the pdf for the power- 
based work estimate p'^j (w) with the TEMA work pdf, p^), in the presence of (. N ) 
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measurements of power during the protocol. Note that while computing p^J, we sum 
over all possible results of the intermediate power measurements. In figures [3] (low initial 
temperature) and [f] (high initial temperature), the cumulative probability for the TEMA 
work pdf with varying number N intermediate power measurements are displayed for 
the LZ problem. For comparison the cumulative probability for the power-based work 
estimate with N = 10 4 measurements is also plotted. In a qualitative way, the two 
distributions approach each other in the limit of N —> oo, but the KS distance (shown in 
the insets) increases with N. This apparent contradiction can be resolved by comparing 
the N —> oo limit of the TEMA work pdf with intermediate power measurements for 
the LZ problem given by 


P\e(w) = 2 


1 + 


1 - 


VT 


2 E„ 


VT 

2 


[PgS(w ~ E 

max ) +p e 5(w + E 

max )] (29) 


S(w) 


with the equivalent expression for the power-based estimate (27). In the low temperature 
case (hence p g ~ 1) depicted in figure [3j the KS distance is maximised at large N, 
because the locations at which the cumulative probabilities perform the largest jumps 
differ. For the power-based estimate, from (27), the jump occurs at w = ur/2 whereas 
for the TEMA based work the jump occurs at w — E max = 2a/( ur/4) 2 + A 2 . Hence 
even for N —y oo, only in the limit that A/(ur) —> 0, the power-based work and the one 
from TEMA with intermediate measurements agree. Note that in the limit A/(ur) —» 0, 
the power operator always commutes with the hamiltonian. 


3. Joint statistics of work from TEMA and power measurements 


To further elucidate the differences between the TEMA and power-based work estimates, 
we consider a modified set-up that allows the simultaneous study of both approaches. 
In order to combine these two approaches, we consider a thought experiment where in 
addition to N power-measurements of the type described in the previous section, also 
energy measurements at the beginning and the end of the force protocol are performed 
according to the TEMA scheme. The outcome of energies e n (0) and e m (r) at the 
beginning and the end of the force protocol and of a sequence x = (x ai , x a2 , • • •, x aN ) of 
eigenvalues of the generalized coordinate X at the equally spaced times of measurement 
H, ,■■■ An occurs with the joint probability (m, x, n ) given by 


^ A (m,x,n) = Tr n m (r)I4(x)n n (0)p 0 n„(0)l/( x ) 

= \{m;r\ U N+1 \Va N )\ 2 |(<A*i| U x |n;0)| 2 

N—l 

n I(</W| U k+1 I<aJ| 2 (n;0| po |ra;0), 

k=l 


X 


(30) 
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Figure 5. Correction factor (e /3 ^ e w p'))^ in the integral fluctuation theorem (|4l|) 


as a function of the number N of power measurements for the LZ problem. The red 
dotted curve is for the parameters v = 5A 2 /h and the black dashed line is for a faster 
sweep rate of v = 40A 2 /h with r = 20A/H. In both cases the initial temperature is 
low and satisfies f3E max = 10. The solid red (thick) and black (thin) lines represent the 
correction factor computed in the N —> oo limit for the slow and fast sweep respectively. 



Figure 6. Correlation between the TEMA work w e and power-based estimate w p for 
the LZ system is displayed as a function of the sweep velocity vt/ A. The sweep time 
is fixed at t = 20H/A. The initial temperature, /3-E max (ur/A) = 10(low temperature, 
red solid) and /3£' max (ur/A) = 0.1 (high temperature, black dashed line), is scaled to 
ensure the same population distribution at different v. 


where, as before, II n (f) denotes the projector on the eigenstate |n;f)[^]of H[X(t)} 
belonging to the eigenenergy e n (f), and V\ is defined in ( flTj ). 

This joint probability and the according probability for the time-reversed process, 
^^(m,x, n), satisfy a detailed balance-like relation of the form 

£P\(jn,x.,n) — e~ /S( ' AF+en ^~ ern ^' > x, m), (31) 

* We do not allow for degeneracy of the energy eigenvalues for the sake of notational simplicity. 
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where x = (xn,Xn- i, ■ ■ •, X\). It holds under two conditions. First, the Hamiltonian H 0 
and the generalized coordinate X must be time reversal invariant meaning that both 
operators commute with the anti-unitary time reversal operator 0, [Hq, 0] = [A", 0] = 
01 This property allows one to express the time evolution running in the backward time 


direction in terms of a propagator proceeding along the forward arrow of time 15,47 
in the following way: 

E^(A) = B t ^r- S ,r-t(A)0, (32) 

The second condition requires that the forward and the backward processes start from 
canonical equilibrium states o) and e~^ H ^ T ^/Z(r), respectively. For a 

detailed derivation we refer to [501. In terms of the conditional probability 

x, n)e^ en ^~ F ^ (33) 

with the free energy F( 0) of the initial state and the analogous conditional probability 
for the backward process 

= &x(n,5c,m)e^ em ^~ F ^\ (34) 

one obtains the microcanonical, detailed balance-like relation[ffl 

x|n) = x|m) (35) 

generalizing the detailed balance-like relation that holds for the forward and backward 
conditional probabilities of energy in the absence of intermediate measurements [48]. 
The joint pdf V^(w e ,w p ) for the two work estimates immediately follows as: 

N 


T’a (w e , w p ) = ^ S ( w e - (e m (r) - e n (0))) 5 [ w p - ^ X(ti)x ai h 

1=1 


n,m,{ai} 

x (m, x, n). 

With the according expression for the time-reversed process, 


N 

V*(w e ,w p )= ^2 0 ■ !r - + ( e m( T ) ~ e„(0))) « ( - V A (t t )x ai h 


(36) 


n,m,{oti} 

x (n, x, m), 


i =1 


(37) 


in combination with the detailed balance-like relation (31) we find a Crooks-type 
fluctuation theorem for the joint distribution of TEMA and power-based work: 

Pa ( -We, ~W p ) = e -h(.w R -AF)j)N (38) 

In full agreement with the earlier observation that any number of intermediate projective 
measurements leaves the Crooks relation, and, consequently, the Jarzynski equality, 

U This condition can be modified for Hamiltonians depending on fields changing their parity under 
time reversal. For details see e.g. 47 48 

ff If the energy eigenvalues e n (t) are degenerate the respective multiplicities d n (t) have to be taken into 
account as ^^(m, x|n)d n (0) = 


48 
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unchanged while modifying the work statistics 49, 50 . One finds for the marginal 


distribution = f dw p Vj?(w e ,w p ) for the TEMA work 

Oe) = e-« AF -»->P0- me ). (39) 

i.e. the Crooks relation being fulfilled. Note that the presence of the unnoticed power 


measurements do have their impact rendering p\ e {w) as defined in (22) different from 


Vfl e (w e ). Yet both work distributions satisfy the Crooks relation in agreement with 

mm- 

For the marginal power-based work pdf Vfl p (w p ) = f dw e Vfl(w e , w p ) only modified 


fluctuation theorems can be obtained. Integrating (38) over all values of w e one obtains 
e -h(AF- Wp ) v N^_ Wp) = V^ p (w p )(e-^- w ^\wX, (40) 

where (-\w p )^ = f dw e ■ Vj?( w e ,w p )/7 D \ tP (w p ) denotes a conditional average. This 
conditional average of the exponentiated work difference determines the modihcation 
relative to the strict Crooks relation. In general it depends on the force protocol. 

Changing A into A and performing an integration over all power-based work values, 
one obtains an integral fluctuation theorem in the form of a modified Jarzynski equality 
reading 


~/3w p \N _ -0AF/ -0(w e -w p )\N 


(e- pw *) N K = e~ p (■ 


A > 


(41) 


where ( • )a = I dw e dw p ■ (w e ,w p ) denotes the average over the joint pdf with N 
power measurements. The correction factor is now determined by the full average of the 
exponentiated difference between TEMA and power-based work. In general, it is also a 
protocol dependent factor. In the limit of a large number of power measurements, we 
can use the approach in section [2] and in a straightforward manner show 

Tr e -/3[A(r)-A(0)]X e -/3H[A(0)] 

' A = W) ‘ 


( e -/30e-™ P ))00 _ 


(42) 


Hence as we remarked in (19) if [H 0 , X ] = 0 is satisfied, the two estimates of work 
w e and w p coincide in the limit of a large number of power measurements (Zeno limit) 
and the correction factor is unity. Figure [5] depicts this correction factor for the LZ 
problem for two speeds of sweeping (see figure captions for details). The correction 
factor behaves non-monotonically with N and can be greater or less than 1 for small N. 
For large N it tends to the correct limiting value as shown in the figure. Moreover for 
the diabatic sweep v = 40A /h, the correction factor is in general significantly smaller 
than the moderate sweep rate case of v — 5A /h. This can be anticipated since in the 
diabatic sweep case with large v, except for a small interval around t — 0 the total 


LZ Hamiltonian (24) approximately commutes with the power operator as the diagonal 


terms ±vt /2 dominate the off-diagonal coupling A. 

So far in this section we elucidated some detailed and integral fluctuation 
relations satisfied by joint TEMA and power-based work estimates and their respective 
marginals and also pointed to formal differences between these two estimates. One 
quantitative measure of the relation between the two estimates of work is the mutual 
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information A I 


A[W e 


W r, 


= In 


( We,Wp) 

'PFM'PjtJvp)' 


While this measure locally quantifies 


A,e\ e /' A 

the interdependence of w e and w p , we consider, the correlation function, as a global 
measure. It is defined as: 

aWe,Wp , (43) 

^"7 /!„ 


Cor(u; e , w p ) = 


where the covariance and standard deviation of the two variable are given by 


a 


We j'tVp 


crz 


= {w e W p )2 ~ {w e )A(w p }%, 

= ( w j)a ~ (( w j)a) 2 ;j=P,e- 


Here, the averages are performed over the joint pdf given by (36). The formal expressions 


for the correlation in general are not very transparent and instead it is more illuminating 
to consider a specific example. We choose the LZ problem introduced in section [2] and 
fix the number of power measurements at N = 100. With a fixed time of sweep at 
r = 20H/A, we ask how the correlation between the two estimates of work varies as the 
velocity of sweep v is changed. In figure [6j we present our results for two thermal initial 
states. For each sweep rate the temperature is chosen such that f3E max (v, A) = {10, 0.1} 
corresponding to fixed populations of the ground state independent of the sweep rate. 
In figure [6] we can clearly see that as vt/ A is increased, the two estimates become more 
and more correlated. This can again be understood as the effect of making the power 
term ycq much larger in magnitude compared to the time-independent part A <j x . This 
effectively renders the commutator between [H 0 ,X] ~ 0 for most of the interval and 
hence the power-based work estimate agrees well with the TEMA based one. Secondly 
we also see that at the smaller value of /3E max corresponding to a larger temperature, 
since the system is more “classical”, the correlation between the two estimates is better. 

In the next section we consider weak continuous measurements of power to estimate 
the work as opposed to projective measurements considered thus far in the paper. 


4. Weak continuous measurement of power 

In the previous sections we modelled the monitoring of the supplied power by means of 
projective measurements of the generalised coordinate X. As we have shown, in the limit 
of large number of such measurements, the system dynamics is frozen in the basis of the 
generalized coordinate and the unitary dynamics generated by the driving plays no role. 
The work statistics in this limit also differs from the ones determined by the TEMA. In 
an attempt to mitigate this situation we consider a weak continuous measurement of the 
generalised coordinate and determine power from such a measurement. One might hope 
that in this case although the estimate of work will be affected by errors inherent in a 
weak measurement process, the measurement backaction will not be so overwhelming 
as to render the unitary dynamics moot. To this end in what follows we will use the 


theory of continuous quantum measurement developed in 52 for our specific situation 
of determining work statistics from power measurements. 
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Before considering a particular scenario of continuous coordinate measurements in 
more detail, we shortly discuss the measurement of the relevant coordinate by means 


of Gaussian Kraus operators 53 as a particular example of a weak instantaneous 
measurement. More general Kraus operators generating a positive operator valued 
measure (POVM) of the coordinate are discussed in the appendix |A} 


4-1. Instantaneous Gaussian coordinate measurements 

We now assume that the measurement of the generalized coordinate X is performed 
with measurement operators consisting of weighted sums of eigen-projection operators 
rather than of a single one. Choosing Gaussian weights we have 


M n = 


(27nr 2 
1 


,V 4 


exp 


\a 


Xf 


4a 2 


(44) 


(27TCT 2 


\ 1/4 


E 


exp 


a 


X r , 


4cr 2 


n 


X 


where a denotes the pointer state indicating the measured value of the coordinate, (Note 
that this choice is slightly different from the one adopted in Appendix [A] where we 
assume that the set of pointer states consists only of the eigenvalues of the generalized 


coordinate in contrast to the continuous range of a values in (44)) and cr 2 the variance 
of the error distribution of measured coordinate values. The work that can be 
estimated from N generalised measurements of X equally spaced in time, is given by 


w, 


= hX ( t n ) a n and its pdf takes the form (see appendix [a] for details) 


N 


(N), ^ 

Pa p ( w ) = 


da k <5 


w — E hX (4) a k 


Tr M ^Mpo, 


(45) 


k =1 


where At = M ajv (tjv)M ajv _ 1 (tjv_i) ■ ■ ■ M ai (ti). In the above, we have used the 
Heisenberg picture representation of the Gaussian measurement operators introduced 
in the appendix (A.5). A natural way to define continuous measurements in this 


framework is to consider weak measurements characterized by a variance that increases 
proportionally to the inverse of the time h between two measurements, i.e. as cr 2 = 
1/(8 kK) where n > 0 quantifies the measurement strength. Hence, with an increasing 
number of measurements less and less information is gained from a single measurement. 


When we take such a scaling for the variance, the measurement operator (44) can be 
expanded as 


M,, 


oc exp 


-2 nh(a — Xfl 


= I 


2 nh(a- X) 2 + 0(h 2 ). (46) 

In appendix [A] we take the generalised measurement operators to be independent of 
the time step. Since the unitary operators can always be written as U k = I + 0(h), 


to lowest order in h we obtain an expression for the pdf (A.6) that depends only on 


the measurement operators and the initial state of the system, provided that the time- 
independent part of the Hamiltonian, Hq, is symmetric in the X basis. In contrast, 
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for the weak measurement case (46) implies that both the measurement and unitary 
operators enter (45) at the same order in h and the system’s state is unchanged to the 
lowest order. In this situation to the best of our attempts, we are not able to determine 
a simple expression for the pdf (45) in the h —> 0 limit. In the next subsection we 
follow another path to tackle this problem by focusing on the differential equation that 
describes a system that is under continuous measurement and is simultaneously driven 
by a change in the Hamiltonian. For the sake of completeness we note that although we 
cannot find a simple closed form expression for the pdf (45) in general, for the trivial 
case of the system Hamiltonian commuting with the measured coordinate [H 0 ,X] = 0, 
we can easily compute the pdf as 

An 

PXp' 


(oo) / \ 

^ J [w) = 


irf 0 T dtA(tfl 


Tr 


exp 


An 


Jo dt\(t)‘ 


w — [A(r) — A(0)] X 


Po 


(47) 


The only difference to the power-based work pdf for projective measurements, 0 is the 
replacement of the delta-function by a Gaussian weight under the trace with a variance 
depending on the force protocol. The average of the exponentiated work yields the 
expression 

(e-P**) = e^^^ 2 (b Tre -/3[A(r)-A(0 )]X po = e ^j; *A*(t) e -/9AF ( 4 g) 


where the second equality on the the right-hand side follows from the commutation of 
Ho and X. Thus we find that even for the trivial situation with [Hq,X] = 0 work 
defined by the integral of weak continuous measurement of power does not satisfy the 
Jarzynski relation. 


4-2. Continuous measurements 


In a final attempt we shall relax the assumption inherent in both the projective and 
the generalized measurement approach that a measurement consists in an instantaneous 
event that interrupts the unitary dynamics of the system. For that purpose we adopt a 
model proposed in 35 . The time evolution of the density matrix caused by the unitary 


dynamics in the presence of continually performed measurements of an observable X is 
described by the following non-linear stochastic master equation (SME) 


m = - j: (ff[A(f)],p(f)] - k[X, (X,p(f)]] (49) 

+ (xpw+ P (t)x - 2 (x(tMo) m, 

where £(£) denotes Gaussian white noise of unit intensity, i.e. (£(£)£(s)) = S(t—s). Apart 
from the first, Hamiltonian term on the right-hand side, the second and third terms are 
supposed to describe the impact of continuous measurements of the observable on the 
system’s dynamics. The second term models the average influence of the measurements, 
and the third, nonlinear and random term accounts for the influence of an individual run 
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Figure 7. Cumulative probability estimated from continuous measurement of power 
in the LZ problem with v = 5A 2 /h, r = 20S./A (blue solid line) and measurement 
strength n = 0.001A/A The initial temperature is small /3E max — 10. The results 
were obtained from 10000-trajectory simulations of the SME (491. The shaded region 
represents the error in the computed pdf represented by the solid line (see text). The 
dashed blue line represents the pdf computed by TEMA. 



w/A 


Figure 8. Cumulative probability of work estimated from continuous measurement of 
power in the LZ problem with a large measurement strength k = 4A /h. The system 
parameters are chosen as v = 0.01A 2 /fi,, r = 20S./A to correspond to an adiabatic LZ 
sweep. The initial temperature is large f3E max = 0.1 for the red curve and smaller 
/3F/ ma x = 10 for the blue curve. 


of measurements. The meanfield-like non-linearity 2 (X)p guaranties the conservation 
of the normalization of the density matrix in the presence of the measurement-induced 
fluctuations. The average ( X ) is performed with respect to the fluctuating density 
matrix p and therefore it is itself a random quantity. In appendix |Bj for the sake of 
pedagogy, we reproduce a derivation of (49) following [35]. There we utilise Gaussian 


measurement operators with variances that are scaled inversely with the time step like 
earlier in this section. 
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Figure 9. Cumulative probability of work estimated from continuous measurement of 
power in the LZ problem with a large measurement strength k = 4A/H. The system 
parameters are chosen as v = 40A 2 /A r = 20H/A to correspond to a diabatic LZ 
sweep. The initial temperature is large /3E max = 0.1 for the red curve and smaller 
/3-Emax = 10 for the blue curve. 


According to the theory of continuous measurements 35 , the result a(t ) of the 
measurement of X at the instant of time t is given by 


a{t) = (X(t)) + -/=£(£), 

V8k 

and consequently, the work based on the measurement of the power becomes 


(50) 


w, 


= / dt\(t)a(t) 


m (<*>« + 


sCi dt. 


(51) 


The SME is understood in the Ito sense. Therefore, the third term on the right-hand 


side of (49) disappears upon an average over realizations of the Gaussian white noise 


yielding for the averaged density matrix p the linear master equation 

§ = -jPWOl.pl-dTh.p]]. 


(52) 


It describes the state of the system in a non-selective measurement of the coordinate X. 
Before discussing the non-linear SME by means of numerical simulations of an 


example we shortly mention a linear SME that is obtained from (49) by disregarding the 
mean-field-type contribution proportional to the average (X) 54,55 . As a consequence, 


the normalization of the density matrix is no longer conserved. Even this linear equation 
is rather complicated for our scenario due to the explicit time dependence of the 
Hamiltonian. In general we find that it can be solved analytically only if Hq and X 
commute. The resulting work pdf coincides as expected with that of the continuous 


action of weak instantaneous Gaussian measurements given by (47). 
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As an example, we generated 10000 realizations of the density matrix by solving 
the SME in (49), displayed in figures [7| [8] and [9] for the LZ model (24) in the presence of 
continuous measurements of a z . For this purpose we used an implicit stochastic Runge- 


Kutta scheme of order 3/2 (56,57 and checked for numerical convergence of single 
trajectories using the method of consistent Brownian paths 57 . We also checked that, 


for the chosen time step, the density matrix remains normalised to 1 to high accuracy 
at all times. We obtained the work pdf by computing a total of 10000 trajectories. For 
the LZ problem the work follows from (51) for each trajectory as 

T 

dt{a z )(t) + 


Wp = 2 


iyo 


V / 8k , 


dt£(t) 


(53) 


Evidently, the work values in a continuous measurement of cr z span the entire real line 
due to the additive contribution of the integral of Gaussian white noise. In order to 
estimate the work pdf from a finite number of simulations, we introduced a binning of 
the size A w p = vh with the same time-step h as used in the discretisation of the SME 
leading to the same set of possible work values as in (26) for A = r/h—1 measurements. 
Based on the independence of different trajectory simulations we could estimate the 
inherent statistical error of the work histograms as the sample variance of 100 blocks of 
100 trajectories each. 

When the measurement strength k is small such that kt <g. 1, the work estimate 
(53) is dominated by the noise term and the resulting pdf is approximately Gaussian 
distributed with zero mean. This feature is independent of the other system parameters 
such as v and A. In figure [7] we plot the cumulative probability for the power-based 
work computed from solving the SME for the same LZ parameters as in section [2j 

When kt > 1, in the work expression the contribution from ( a z )(t ) dominates 
over the noise. Thus the behaviour of the cumulative probability in figures [8] and [9] 
can be well understood from the behaviour of ( a z )(t ). For instance, from (53) this 
immediately explains why the range of allowed work values is comparable to the interval 
[—nr/2, nr/2]. Also note that for our choice of parameters in figures [8] and [9] we are also 
in the strong measurement regime of k > A/h for the LZ problem where the coherent 
dynamics rate A is trumped by the measurement backaction. It is known 58-60 that 
for the LZ system under continuous measurement, the behaviour of single trajectory 
solutions of (49) goes from near unitary at very small measurement strength to the 
so-called “random-telegraph” dynamics at the strong measurement k > A/h regime 
(In (58,601 the dynamics in the strong measurement regime is referred to as the quantum 
Zeno effect. In our work we reserve the latter term for the complete freezing of unitary 
dynamics achieved by repeated projective measurements discussed in section |2]). The 
random-telegraph behaviour is characterized by the population difference of diabatic 
states (a z ) remaining localised either at ±1 and undergoing rapid transitions between 
the two values at different times during the evolution (see discussion in Sec. Ill B 


of [58]). In the adiabatic sweep case depicted in figure |8j since vr <C A the initial 
state is off-diagonal (with almost equal distribution amongst the ±1 eigenstates) in 
the diabatic i.e. a z basis and throughout the dynamics the system Hamiltonian has 
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a significant off-diagonal component in the diabatic basis i.e. [H[X{t)], a z \ f 0. Since 
there is equal contribution from both (a z ) = ±1 to (53), the work pdf is centered 
about w — 0 in figure [8j Nonetheless the power-based cumulative probability does not 
capture the sharp jump at w = 0 of the TEMA estimate. In the diabatic sweep in figure 
[9j the system Hamiltonian and the initial state are almost diagonal in the diabatic 
basis giving cumulative probabilities that have sharp jumps at w = nr/2 (low initial 
temperature) or w = ±ur/2 (high temperatures). Also in figure [9] the jumps in the 
power-based and TEMA estimates occur at approximately the same values of w. We 
note that this behaviour owes its explanation, as noted in section [2j to the fact that for 
our choice vr /4 3> A the magnitude of maximum work possible in TEMA E max vr/2 
almost coincides with the location of the jumps in the power measurement cumulative 
probability. 

In summary we find that the work pdf estimated even by weak continuous 
measurement of the power does not reproduce the behaviour of the TEMA estimate in 
general. For very weak measurement strengths nr <C 1, the estimate is noise dominated 
and has Gaussian behaviour irrespective of v and A. When kt > 1 (strong measurement 
regime) the pdf depends on the sweep rate v (for a fixed value of r) which sets the extent 
to which the total Hamiltonian is off-diagonal in the measurement basis. 


5. Conclusion 

This work serves as a detailed illustration of the difficulties involved in defining work 
in quantum systems in a manner analogous to classical systems. We considered the 
statistical properties of work, defined as the integral of supplied power, for a quantum 
system. We showed that even a careful definition of work in terms of repeated 
measurements of the system’s instantaneous power leads to a statistics quite different 
from the usual two energy measurement approach (TEMA) for defining work. In the 
power-based approach in general even the Jarzynski equality or the Crooks relation 
do not hold. In the limit of a large number of projective measurements of the 
instantaneous power we found that the quantum Zeno effect leads to a freezing of the 
system’s dynamics in the power operator’s basis. The statistics of power-based work 
in this limit is very different from TEMA except for the trivial case when the power 
operator commutes with the system Hamiltonian at all times. Furthermore we carried 
out a detailed comparison of the power-based work and the TEMA work estimate by 
considering the joint pdf of both types of work estimates in a setting that combines both 
approaches. We obtained a Crooks type fluctuation theorem for the joint probability 
distribution and modified integral fluctuation theorems for the marginal power-based 
work distribution. We also studied the correlation between the two estimates of work 
using the joint probability for the LZ problem. Finally, relaxing the condition of 
projective, instantaneous measurements, weak continuous measurement of power was 
discussed within the formalism of stochastic master equation (SME). Using the Landau- 
Zener problem as an example, we determined the pdf of power-based work numerically 
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and analysed its properties in the limiting cases of weak and strong measurement 
strengths. We found that since the power operator in general does not commute with the 
total Hamiltonian of the system, also within the framework of continuous measurements 
the power-based estimate is not able to reproduce the qualitative features of the TEMA 
work estimate. 

Finally we would like to note that the treatment of continuous measurements in 
terms of an SME also leads to a Lindblad-type master equation, see (52), similarly as 
obtained in various recent publications 61-64 which are concerned with the definition of 
work and heat and the quest for fluctuation theorems in open quantum systems. There, 
the presence of terms in the equation of motion of the reduced density matrix describing 
energy-non-conserving effects is caused by the interaction with the environment. In 
contrast, in the present case it is solely the result of a weak but continuously acting 
measurement that leads to similar formal structures. 
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A. Generalized power measurements 

In this appendix we consider generalized measurements of the power observable during 
the protocol and prove a bound on the exponential average of the work estimate obtained 
in this manner. We replace the projective measurements of the coordinate X in the 
definition of joint probability ( [TT| by a set of self-adjoint, i.e. minimally disturbing 
generalized measurements which have the following form 

= Y y/pfaM n i- (A- 1 ) 

Oik 

We associate with each outcome a 3 of the generalized measurement the coordinate 
eigenvalue x Qj . p^aflaQ, then denotes the conditional probability for erroneous 
assignment of eigenvalue x aj to x ak by the generalized measurement. These must be 
real and positive and in accordance with partition of unity, ^2 aj Aft. M a . = /, they add 
up to unity: 

Y'P( a j\ a k) = 1 . 


(A.2) 
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Thus, in the same manner as in section [2] we can write down the pdf for the work 
estimated by summing up results of generalized measurements of power as: 

( N 

P [ Kp\ W ) = X 5 Iw-^X^X^h 

K} V i= 1 


j Tr VKM (x) PoVl M (x), (A.3) 


where 


V\m ( x ) — U} Njto M olN UNM olN _ 1 ■ ■ ■ U2M ai Ui 
— M olN (t N )M aN _ 1 {t N -i) ■ ■ ■ M ai (ti), 


(A.4) 

(A.5) 


where in the second line we have defined Heisenberg picture measurement operators 
M a .{tj) = Ul. to M a .Ut jt t 0 - 

In the limit of large N, in line with the treatment in section [2j the unitary operators 
) can be expanded as Uj — I + 0(h). Provided the time-independent part of 
the Hamiltonian H 0 is symmetric in the coordinate X’s basis i.e. {p aj \ H 0 \ip ak ) = 
(<fa k | Ho | (p aj ), the individual terms that are of order h cancel amongst themselves. 
This ensures that such terms do not sum up to a term that adds to the leading order 
contribution. Hence to leading order in h the work distribution is given by: 

( N 

(ad V <=1 

+ 0(h), (A.6) 


TrMl.Mi. . ■■■Ml o 0 


a N 


ajv-l 


in (A.4 


where we have used the fact that generalized measurement operators at different times 
commute. In the trivial case when the Hamiltonians at different times commute, the 


unitary operators Uk commute with the measurement operators and (A.6) becomes 
exactly valid not just when h —» 0. 

Unlike for the case of projective measurements of X, further simplification of the 


expression (A.6) by taking the limit A^ —)• oo is not possible without choosing a specific 


form for the error distribution function p(aj\ak). Nonetheless as we show below, we 
can derive an inequality for the exponential average under some assumptions. The 
exponential average of the work estimate to zeroth order in h is given by: 

N 

(e- 0w )M = Y, M Po M , (A.7) 

k i= 1 


where ( f(X))k,M = J2 aj f( x aj)p( a j\ a k)- Applying the Jensen inequality with respect 
to the measurement distribution, ( e A )k,M > e^ k ' M , to each term in the product on the 
right-hand side of ( |A.7[ ) and taking the limit h 0 we get 

(e^ 0w ) M > 5>P (~/3(A (r) - A(0)) (<p ak \po \Va k ) ■ 

k 

Finally if we assume that we have a homogeneous distribution for p(aj\aQ, dependent 
only on the difference of the eigenvalues {x aj , x ak } and centered about x ak for each k, 
the mean (X) k ^ M will equal x ak leading to 

(e~P w ) M > X e xp(-/3(A(r) - A(0))x a J (<p a J p 0 \<p ak ) ■ 

k 
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Comparing the above equation to (18), we have the inequality 


(A.8) 


that is the exponential average of power-based work obtained from a generalised 
measurement of power is always greater than the one obtained from projective power 
measurements. 


B. Continuous measurement stochastic master equation 


Let us consider the state of a system p(t) that is normalized at time t and is subject 


to the Gaussian measurement introduced in (44) and with outcome cq and a unitary 


evolution operator lJ t for a short interval h. We also chose the variance of the Gaussian 
measurement as a 2 = 1/(8 kK). The normalised state at t + h is given by 

M at U t p(t)U}M at 


p(t + h) — 


Tr (M at U t p(t)U} 


(B.l) 


atj 


Examining the probability distribution for the measurement outcome, given by the 


denominator in (B.l), as shown in 35 one can write (provided h is chosen small enough 
so that the Gaussian is broader than the X distribution of the state p(t)) 

AW 

w (R2) 

where AW is a zero mean Gaussian random variable with variance h and ( X ) = 
Tr (p(t)X). Taking into account the scaling AW 2 ~ h and (B.2), we expand the 


operators in (B.l) as 


U t =I-i-H t + 0(h 2 ), 


(B.3) 


M n 


cX / - k(2 h - AW 2 ) (X - (X)) + V2kAW (x - (X)) 

+ 0(h 3/2 ), (B.4) 


where we have denoted the instantaneous Hamiltonian as H t . Substituting (B.3) and 


(B.4) we get 


h 


p(t + h) = pit) - i-[H t ,p(t)] + V2fAW |X - (X), 


(B.5) 


k(2 h - AW 2 ) jx 2 , p(t)j + 2kAW 2 Xp(t)X 
Ak{AW 2 - h) Tr (X - (X)Y p{t) + 0{h 3/2 ). 


Taking the infinitesimal limit h —y dt and AW —y dW = £dt, the Wiener differential in 
(B.5) and observing that dW 2 = dt, we immediately obtain (49). 
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